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Abstract 

Cosmological inflation in supergravity is a great window to High-Energy Physics be- 
yond the Standard Model. I review a construction of the new F(1Z) supergravity the- 
ories, and consider some of their applications to cosmology. The F(JZ) supergravity 
theories are the N = 1 locally supersymmetric extensions of the well studied f(R) 
gravity theories in four space-time dimensions, which are often used for 'explaining' 
inflation and Dark Energy. A manifeslty supersymmetric description of F(TZ) super- 
gravities exists in terms of N — 1 superfields, by using the (old) minimal Poincare 

(N 



supergravity in curved superspace. We find that any F(1Z) supergravity is classically 
equivalent to the Poincare-type matter-coupled N = 1 supergravity via the super- 
O ■ field Legendre-Weyl-Kahler transformation. The (non-trivial) Kahler potential and the 

scalar superpotential of the inflaton superfield are determined in terms of the original 
holomorphic F(JZ) function. The conditions for stability, the absence of ghosts and 
tachyons are also found. No-scale F{1Z) supergravity is constructed too. Three dif- 
ferent examples of F(TZ) supergravity theories are studied in detail. The first example 
is devoted to recovery of the standard (pure) N = 1 supergravity with a negative cos- 
mological constant from F(1Z) supergravity. As the second example, a generic 1Z 2 
supergravity is investigated, and the existence of the AdS bound on the scalar curva- 
ture is discovered. As the third example, a simple viable realization of chaotic inflation 
in supergravity is given, via an embedding of the Starobinsky inflationary model into 
the F(1Z) supergravity. It is also found that F(1Z) supergravity can have a positive 
cosmological constant in the low-curvature regime. Our approach does not introduce 
new exotic fields or new interactions, beyond those already present in (super)gravity. 
A nonminimal scalar-curvature coupling in gravity and supergravity, Higgs inflation 
and its correspondence to Starobinsky inflation in f(R) gravity and F(1Z) supergrav- 
ity, respectively, are establlished. In the outlook, reheating and particle production in 
the post-inflationary universe, CP- violation, origin of baryonic asymmetry, lepto- and 
baryo-genesis are briefly discussed. F(1Z) supergravity shows promise for possible 
solutions to those outstanding problems. 



CONTENTS 



1 . Introduction and Motivation 3 

2. Starobinsky Model of Inflation 6 

3. Gravity and Scalar-Tensor Gravity 7 

4. Inflationary Theory and Observations 9 

5. Supergravity and Superspace 13 

6. F(1Z) Supergravity in Superspace 14 

7. No-scale F(R) Supergravity 18 

8. Fields from Superfields in F(1Z) Supergravity 20 

9. Generic 1Z 2 Supergravity, and AdS Bound 21 

10. Chaotic Inflation in F(1Z) Supergravity 26 

11. Cosmological Constant in F(1Z) Supergravity 26 

12. Nonmimimal Scalar-Curvature Coupling in Gravity and Supergravity, 

and Higgs Inflation 34 

13. Reheating and Quantum Particle Production 41 

14. Conclusion 43 

15. Outlook: CP- violation, Baryonic Asymmetry, Lepto- and Baryo-genesis, 
Non-Gaussianity, and Experimental Tests 45 

16. Acknowledgements 47 

17. References 48 



This review is based on the original papers [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11] and 
the invited seminar talks given by the author at the CERN Theory Division in Geneva 
(Switzerland), the DESY Theory Group in Hamburg, the Max-Planck Institutes for 
Physics and Gravitational Physics in Munich and Golm, respectively, Heidelberg Uni- 
versity (Germany), The University of Tokyo, Waseda University, Tokyo Institute of 
Technology, Tokyo University of Science, KEK Theory Group and Yukawa Institute 
of Theoretical Physics in Kyoto (Japan). 



2 



Figure 1: Timeline of the universe according to WMAP [12] 

1 Introduction and Motivation 

A very brief history of our universe is given by the standard Fig. 1 taken from the 
website of the Wilkinson Microwave Anisotropy Probe (WMAP) [12]. We focus on 
the field-theoretical description of the inflationary phase of early universe and its post- 
inflationary dynamics (reheating and particle production) in the context of supergrav- 
ity. To begin with, let us first introduce some basics of inflation. 

Cosmological inflation (a phase of 'rapid' quasi-exponential accelerated expan- 
sion of universe) [13, 14, 15] predicts homogeneity of our universe at large scales, 
its spatial flatness, large size and entropy, and the almost scale-invariant spectrum of 
cosmological perturbations, in good agreement with the WMAP measurements of the 
CMB radiation spectrum [12, 16]. Inflation is also the only known way to generate 
structure formation in the universe via amplifying quantum fluctuations in vacuum. 

However, inflation is just the cosomological paradigm, not a theory! The known 
field-theoretical mechanisms of inflation use a slow-roll scalar field <\> (called inflaton) 
with proper scalar potential V(<f>) [14, 15]. 

The scale of inflation is well beyond the electro-weak scale, ie. is well beyond 
the Standard Model of Elementary Particles! Thus the inflationary stage in the early 
universe is the most powerful High-Energy Physics (HEP) accelerator in Nature (up to 
10 10 TeV). Therefore, inflation is the great and unique window to very HEP! 

The nature of inflaton and the origin of its scalar potential are the big mysteries. 

Throughout the paper the units h = c = 1 and the spacetime signature (+,—,—,—) 
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Figure 2: CMB from WMAP [12]: Mona Lisa of Cosmology 



are used. See ref. [17] for our use of Riemann geometry of a curved spacetime. 

The Cosmic Microwave Background (CMB) radiation picture from the WMAP 
(Fig. 2) is one of the main sources of data about early universe. Deciphering this pic- 
ture in terms of density perturbations, gravity wave polarization, power spectrum and 
its various indices is a formidable task. It also requires the heavy CMB mathematical 
formalism based on General Relativity — see eg., the textbooks [18, 19, 20]. Fortu- 
nately, we do not need that formalism for our purposes, since the relevant indices can 
also be introduced in terms of the inflaton scalar potential (Sec. 4). We assume that in- 
flation did happen. There exist many inflationary models — see eg. the textbook [15] 
for their description and comparison (without supersymmetry). Our aim is a viable 
theoretical description of inflation in the context of supergravity. 

The main Cosmological Principle of a spatially homogeneous and isotropic (1 + 3)- 
dimensional universe (at large scales) gives rise to the FLRW metric 

(1) 

where the function a(t) is known as the scale factor in 'cosmic' (comoving) coordi- 
nates (t,r,9,(p), and k is the FLRW topology index, k = (-1,0, +1). The FLRW 
metric (1) admits the six-dimensional isometry group G that is either 5*0(1, 3), E(3) 
or SO(A), acting on the orbits G/SO(3), with the spatial three-dimensional sections 
H 3 , E 3 or S 3 , respectively. The Weyl tensor of any FLRW metric vanishes, 



r/s 2 



dt 2 - a 2 (t) 



1 — kr' 



(2) 



where //, u, A, p — 0, 1, 2, 3. The early universe inflation (acceleration) means 



where H —a /a is called Hubble function. We take k = for simplicity. The amount 
of inflation (called the e-foldings number) is given by 



Next, a few words about our strategy. It is well recognized by now that one has 
to go beyond the Einstein-Hilbert action for gravity, both from the experimental view- 
point (eg., because of Dark Energy) and from the theoretical viewpoint (eg., because 
of the UV incompleteness of quantized Einstein gravity, and the need of its unification 
with the Standard Model of Elementary Particles). 

In our approach, the origin of inflation is purely geometrical, ie. is closely re- 
lated to space-time and gravity. It can be technically accomplished by taking into 
account the higher-order curvature terms on the left-hand-side of Einstein equations, 
and extending gravity to supergravity. The higher-order curvature terms are supposed 
to appear in the gravitational effective action of Quantum Gravity. Their derivation 
from Superstring Theory may be possible too. The true problem is a selection of those 
high-order curvature terms that are physically relevant or derived from a fundamental 
theory of Quantum Gravity. 

There are many phenomenological models of inflation in the literature, which usu- 
ally employ some new fields and new interactions. It is, therefore, quite reasonable 
and meaningful to search for the minimal inflationary model building, by getting most 
economical and viable inflationary scenario. I am going to use the one proposed the 
long time ago by Starobinsky [21, 22], which does not use new fields (beyond a space- 
time metric) and exploits only gravitational interactions. I also assume that the general 
coordinate invariance in spacetime is fundamental, and it should not be sacrificed. 
Moreover, it should be extended to the more fundamental, local supersymmetry that 
is known to imply the general coordinate invariance. It thus leads us to supergrav- 
ity which, in addition, automatically has several viable candidates for Dark Matter 
particle (see Sec. 15 for more). 

On the theoretical side, the available inflationary models may be also evaluated 
with respect to their "cost", ie. against what one gets from a given model in relation 
to what one puts in! Our approach does not introduce new fields, beyond those already 
present in gravity and supergravity. We also exploit (super)gravity interactions only, 
ie. do not introduce new interactions, in order to describe inflation. 

Before going into details, let me address two common prejudices and objections. 

The higher-order curvature terms are usually expected to be relevant near the space- 
time curvature singularities. It is also quite possible that some higher-derivative grav- 
ity, subject to suitable constraints, could be the effective action to a quantized theory 




(3) 




(4) 
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of gravity, like eg., in String Theory. However, there are also some common doubts 
against the higher-derivative terms, in principle. 

First, it is often argued that all higher-derivative field theories, including the higher- 
derivative gravity theories, have ghosts (i.e. are unphysical), because of Ostrogradski 
theorem (1850) in Classical Mechanics. As a matter of fact, though the presence of 
ghosts is a generic feature of the higher-derivative theories indeed, it is not always 
the case, while many explicit examples are known (Lovelock gravity, Euler densities, 
some f(R) gravity theories, etc.) — see eg., ref. [23] for more details. In our approach, 
the absence of ghosts and tachyons is required, while it is also considered as one of the 
main physical selection criteria for the good higher-derivative field theories. 

Another common objection against the higher-derivative gravity theories is due to 
the fact that all the higher-order curvature terms in the action are to be suppressed by 
the inverse powers of M P1 on dimensional reasons and, therefore, they seem to be 'very 
small and negligible' . Though it is generically true, it does not mean that all the higher- 
order curvature terms are irrelevant at all scales much less than M P1 . For instance, it 
appears that the quadratic curvature terms have dimensionless couplings, while they 
can easily describe the early universe inflation (in the high-curvature regime). A non- 
trivial function of R in the effective gravitational action may also 'explain' Dark En- 
ergy in the present universe (see Sec. 1 1 for more). 

2 Starobinsky model of inflation 

It can be argued that it is the scalar curvature-dependent part of the gravitational ef- 
fective action that is most relevant to the large-scale dynamics H(t). Here are some 
simple arguments. 

In 4 dimensions all the independent quadratic curvature invariants are R^ uXp R^uXp , 
R^R^ and R 2 . However, 

J d 4 x V=g (R^ xp R,ux P - 4i^iV + R 2 ) (5) 

is topological (ie. a total derivative) for any metric, while 

J d 4 x (SR^R^ - R 2 ) (6) 

is also topological for any FLRW metric, because of eq. (2). Hence, the FLRW-relevant 
quadratically-generated gravity action is (87rGjv = 1) 

S = ~J d 4 x^(R-R 2 /M 2 ) (7) 
1 To the best of my knowledge, this proposal was first formulated by A.D. Sakharov in 1967. 



6 



This action is known as the (simplest) Starobinsky model [21, 22]. Its equations of 
motion allow a stable inflationary solution, and it is an attractor! When H 3> M, one 
finds 

H ~ (~6~) ^ cnd "^ (8) 

It is the particular realization of chaotic inflation (ie. with chaotic initial conditions) 
[24], and with a Graceful Exit. 

In the case of a generic gravitational action with the higher-order curvature terms, 
the Weyl dependence can be excluded due to eq. (2) again. A dependence upon the 
Ricci tensor can be also excluded since, otherwise, it would lead to the extra propa- 
gating massless spin-2 degree of freedom (in addition to a metric) described by the 
field dC/dR^. The higher derivatives of the scalar curvature in the gravitational La- 
grangian £ just lead to more propagating scalars [25], so I simply ignore them for 
simplicity in what follows. 



3 f(R) Gravity 



The Starobinsky model (7) is the special case of the f(R) gravity theories [26, 27, 28] 
having the action 

1 



In the absence of extra matter, the gravitational (trace) equation of motion is of the 
fourth order with respect to the time derivative, 

3 d ( ( ^m^_ - /;/ . /; , ; . /; :j 



a 



3 dt\ dt 



where we have used H — ^ and R = —Q(H +2H 2 ). The primes denote the derivatives 
with respect to R, and the dots denote the derivative with respect to t. Static de-Sitter 
solutions correspond to the roots of the equation Rf'(R) = 2f(R) [29]. 

The 00-component of the gravitational equations is of the third order with respect 
to the time derivative, 

3/7^^ - 3(H +H 2 )f{R) - l -f{R) =0 (11) 

The (classical and quantum) stability conditions in f(R) gravity are well known 
[26, 27], and are given by (in our notation) 

f(R) > and f"(R) < (12) 

respectively. The first condition (12) is needed to get a physical (non-ghost) graviton, 
while the second condition (12) is needed to get a physical (non-tachyonic) scalaron 
(see Sec. 9 for more). 
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Any f(R) gravity is known to be classically equivalent to the certain scalar-tensor 
gravity having an (extra) propagating scalar field [30, 31, 32]. The formal equivalence 
can be established via a Legendre-Weyl transform. 

First, the f(R) -gravity action (9) can be rewritten to the form 



2k 2 



d A x y/^g {AR — Z(A)} 



(13) 



where the real scalar (or Lagrange multiplier) A(x) is related to the scalar curvature R 
by the Legendre-like transformation: 



R = Z'(A) and f(R) = RA(R) - Z(A(R)) 



(14) 



with k 2 = 8nG N = Mpj 2 . 

Next, a Weyl transformation of the metric, 



9vu,(x) ->■ exp 

with arbitrary field parameter <j>(x) yields 

2k6(x 



2n(j)(x) 
V6 



(15) 



^gR^- V^g exp 



Therefore, when choosing 



A(K(f>) = exp 



-2k<P(x) 
V6 



(16) 



(17) 



and ignoring a total derivative in the Lagrangian, we can rewrite the action to the form 



d 4 x \J—g 



2^ + IgTd^ 



2k 2 



exp 



4:K(j)(x) 



Z{A{k4>)) 



(18) 



in terms of the physical (and canonically normalized) scalar field <p{x), without any 
higher derivatives and ghosts. As a result, one arrives at the standard action of the 
real dynamical scalar field cj)(x) minimally coupled to Einstein gravity and having the 
scalar potential 



v{4>) 



M 2 X 



exp 



MpiV6 



Z ( exp 



-20 



MpiVQ. 



(19) 



In the context of the inflationary theory, the scalaron (= scalar part of spacetime 
metric) can be identified with inflaton. This inflaton has clear origin, and may also 
be understood as the conformal mode of the metric over Minkowski or (A)dS vacuum. 
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Figure 3: The inflaton scalar potential v(x) = [e y — l) 2 in the Starobinsky model, after 

y ->■ -y 

In the Starobinsky case of f(R) = R — R 2 /M 2 , the inflaton scalar potential reads 

V(y) = V (e-y-l) 2 (20) 
where we have introduced the notation 

y = J--^- and V = -MlM 2 (21) 
y V 3 Mpi 8 P1 

It is worth noticing here the appearance of the inflaton vacuum energy V driving in- 
flation. The end of inflation (Graceful Exit) is also clear: the scalar potential (20) has a 
very flat (slow-roll) 'plateau', ending with a 'waterfall' towards the minimum (Fig. 3). 

It is worth emphasizing that the inflaton (scalaron) scalar potential (20) is derived 
here by merely assuming the existence of the R 2 term in the gravitational action. The 
Newton (weak gravity) limit is not applicable to early universe (including its inflation- 
ary stage), so that the dimensionless coefficient in front of the R 2 term does not have to 
be very small at early time. It distinguishes the primordial 'Dark Energy' driving infla- 
tion in the early Universe from the 'Dark Energy' responsible for the present universe 
acceleration. 

4 Inflationary Theory and Observations 

The slow-roll inflation parameters are defined by 

e(0 = ±JH*n^ and r/(0) = M 2 ^ (22) 

A necessary condition for the slow-roll approximation is the smallness of the inflation 
parameters 

£(</>)<! and \r]((j))\ <1 (23) 
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The first condition implies a (t) > 0. The second one guarantees that inflation lasts 
long enough, via domination of the friction term in the inflaton equation of motion, 

3H 0= -V. 

The CMB temperature fluctuations (Fig. 2) have the scale 5T/T ps 10~ 5 at the 
WMAP normalization of 500 Mpc. As is well known [15], scalar (p s ) and tensor (p t ) 
perturbations of metric do decouple. The scalar perturbations couple to the density of 
matter and radiation, so they are responsible for the inhomogeneities and anisotropics 
in the universe. The tensor perturbations (or gravity waves) also contribute to the 
CMB, while their experimental detection would tell us much more about inflation. 
The CMB raditation is expected to be polarized due to Compton scattering at the time 
of decoupling [33, 34]. 

The primordial spectrum is proportional to /c n ~\ in terms of the comoving wave 
number k and the spectral index n, in the 2-point function (observable) 

^mmyjp^^r (24) 

In theory, the slope n s of the scalar power spectrum, associated with the density 

perturbations, ^ fc™ 8-1 , is given by n s = 1 + 2i] — 6e, the slope of the tensor 

primordial spectrum, associated with gravitational waves, is n t = —2e, and the tensor- 
to-scalar ratio is r = Sp s /5p t = 16e (see eg., ref. [15]). 

It is straightforward to calculate those indices in any inflationary model with a 
given inflaton scalar potential. In the case of the Starobinsky model and its scalar 
potential (20), one finds [35, 36, 6] 

2 3\nN e 2 /ln 2 iV e \ 

and 

12 

r w — « 0.004 (26) 

iV e 2 

with N e ~ 55. The very small value of r is the sharp prediction of the Starobinsky 
inflationary model towards r -measurements in a future. 

Those theoretical values are to be compared to the observed values of the CMB 
radiation due to the WMAP satellite mission — see a picture of the WMAP satellite in 
Fig. 4. For instance, the most recent WMAP7 observations [12] yield 

n s = 0.963 ± 0.012 and r < 0.24 (27) 

with the 95 % level of confidence. 

The amplitude of the initial perturbations, A|, = M^V/ (247r 2 e), is also a physical 
observable, whose experimental value is known due to another Cosmic Background 
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Figure 4: The WMAP satellite mission [12] 



Explorer (COBE) satellite mission — see ref. [37] and Fig. 5 for a picture of the COBE 
satellite; 

1/4 

= 0.027 Mpi = 6.6 x 10 16 GeV (28) 
It determines the normalization of the i? 2 -term in the action (7) 

M l~2 p~ v 

— = 4-J-- (2.7) 2 • — — - ■ 10" 4 « (3.5 ± 1.2) ■ lO" 6 (29) 

Mpi V 3 V ; (1 - e-v ) 2 V ; 

In particular, the inflaton mass is given by M inf = MJ a/6. 

The main theoretical lessons, that we can draw from the discussion above towards 
our next goals, are: 

(i) the main discriminants amongst all inflationary models are given by the values of 
n s and r; 

(ii) the Starobinsky model (1980) of chaotic inflation is very simple and economic. It 
uses gravity interactions only. It predicts the origin of inflaton and its scalar potential. 
It is still viable and consistent with all known observations. Inflaton is not charged 





Figure 5: The COBE satellite mission [37] 
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Figure 6: Starobinsky inflation vs. m 2 2 /2 and A</> 4 



(singlet) under the SM gauge group. The Starobinsky inflation has an end (Graceful 
Exit), and gives the simple explanation to the WMAP-observed value of n s . The key 
difference of Starobinsky inflation from the other standard inflationary models (having 
\m 2 (f) 2 or A</> 4 scalar potentials) is the very low value of r — see the standard Fig. 6 for 
a comparison and ref. [38] for details. A discovery of primordial gravitational waves 
and precision measurements of the value of r (if r > 0.1) with the accuracy of 0.5% 
may happen due to the ongoing PLANCK satellite mission — see ref. [39] and Fig. 7 
for a picture of the PLANCK satellite; 

(iii) the viable inflationary models, based on f(R) = R + f(R) gravity, turn out to 
be close to the simplest Starobinsky model (over the range of R relevant to inflation), 
with f(R) ~ R 2 A(R) and the slowly varying function A(R) in the sense 




Figure 7: Ongoing Planck satellite mission [39] 
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5 Supergravity and Superspace 



Super symmetry (SUSY) is the symmetry between bosons and fermions. SUSY is the 
natural extension of Poincare symmetry, and is well motivated in HEP beyond the SM. 
Supersymmetry is also needed for consistency of strings. Supergravity (SUGRA) is 
the theory of local supersymmetry that implies general coordinate invariance. In other 
words, considering inflation with supersymmetry necessarily leads to supergravity. As 
a matter of fact, most of studies of superstring- and brane-cosmology are also based 
on their effective description in the 4-dimensional N — 1 supergravity. 

It is not our purpose here to give a detailed account of SUSY and SUGRA, because 
of the existence of several textbooks — see e.g., refs. [40, 41, 42]. In this Section I 
recall only the basic facts about N — 1 supergravity in four spacetime dimensions, 
which are needed for our main purposes. 

A concise and manifestly supersymmetric description of SUGRA is provided by 
Superspace. In this section the natural units c = h = k = M P1 = 1 are used. 

Supergravity needs a curved superspace. However, they are not the same, because 
one has to reduce the field content to the minimal ones corresponding to the off- shell 
supergravity multiplets. It is done by imposing certain constraints on the supertorsion 
tensor in curved superspace [40, 41, 42]. An off-shell supergravity multiplet has some 
extra (auxiliary) fields with noncanonical dimensions, in addition to physical spin-2 
field (metric) and spin-3/2 field (gravitino). It is worth mentioning that imposing the 
off-shell constraints is independent upon writing a supergravity action. 

One may work either in a full (curved) superspace or in a chiral one. There are 
some anvantages of using the chiral superspace, because it helps us to keep the auxil- 
iary fields unphysical (i.e. nonpropagating). 

The chiral superspace density (in the supersymmetric gauge-fixed form) reads 

£(x, 9) = e(x) [1 - 2i6a a ^ a (x) + 9 2 B(x)] , (31) 

where e = det g^ u , g^ u is a spacetime metric, ip^ = e^ip£ is a chiral gravitino, 
B = S — iP is the complex scalar auxiliary field. We use the lower case middle greek 
letters /i, u, . . . = 0, 1, 2, 3 for curved spacetime vector indices, the lower case early 
latin letters a, b, . . . = 0, 1, 2, 3 for flat (target) space vector indices, and the lower case 
early greek letters a, (3, . . . = 1,2 for chiral spinor indices. 

A solution to the superspace Bianchi identities together with the constraints defin- 
ing the N = 1 Poincare-type minimal supergravity theory results in only three covari- 
ant tensor superfields 1Z, Q a and W^, subject to the off-shell relations [40, 41, 42]: 

Ga = Qa, Wafh = W {afh) , V.ft = V. W a p~ f = , (32) 

and 

= V « n > V7W ^ = %Z + , (33) 

where (V , V..V .) stand for the curved superspace N — 1 supercovariant deriva- 
tives, and the bars denote complex conjugation. 
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The covariantly chiral complex scalar superfield 1Z has the scalar curvature R as the 
coefficient at its 9 2 term, the real vector superfield Q , has the traceless Ricci tensor, 
R^iu + R»n — hdnvR, as the coefficient at its 9a a 9 term, whereas the covariantly chiral, 
complex, totally symmetric, fermionic superfield W a/ 3 7 has the self-dual part of the 
Weyl tensor C a/ 3 y s as the coefficient at its linear ^-dependent term. 

A generic Lagrangian representing the supergravitational effective action in (full) 
superspace, reads 

c = c(n,g,w,...) 04) 

where the dots stand for arbitrary supercovariant derivatives of the superfields. 

The Lagrangian (34) it its most general form is, however, unsuitable for physi- 
cal applications, not only because it is too complicated, but just because it generi- 
cally leads to propagating auxiliary fields, which break the balance of the bosonic and 
fermionic degrees of freedom. The important physical condition of keeping the su- 
pergravity auxiliary fields to be truly auxiliary (ie. nonphysical or nonpropagating) 
in field theories with the higher derivatives was dubbed the 'auxiliary freedom ' in 
refs. [43, 44]. To get the supergravity actions with the 'auxiliary freedom', we will use 
a chiral (curved) superspace. 



6 F(1Z) Supergravity in Superspace 

Let us first concentrate on the scalar-curvature- sector of a generic higher-derivative 
supergravity (34), which is most relevant to the FRLW cosmology, by ignoring the 
tensor curvature superfields W Q/ g 7 and Q ., as well as the derivatives of the scalar 
superfield 1Z, like that in Sec. 2. Then we arrive at the chiral superspace action 

S F = J d A xd 2 0£F{TZ) + H.c. (35) 

governed by a chiral or holomorphic function F(1Z). 2 Besides having the manifest 
local N = 1 supersymmetry, the action (35) has the auxiliary freedom since the aux- 
iliary field B does not propagate. It distinguishes the action (35) from other possible 
truncations of eq. (34). The action (35) gives rise to the spacetime torsion generated 
by gravitino, while its bosonic terms have the form 

S f = -\j d 4 x^f(R) (36) 

Hence, eq. (35) can also be considered as the locally N — 1 supersymmetric exten- 
sion of the /(i?)-type gravity (Sec. 3). However, in the context of supergravity, the 
'supersymmetrizable' bosonic functions f(R) are very restrictive (see Sees. 9 and 10). 

2 The similar component field construction, by the use of the 4D, N = 1 superconformal tensor 
calculus, was given in ref. [45]. 



14 



The superfield action (35) is classically equivalent to 

S v = J d 4 xd 2 d£[ZK-V(Z)} +H.c. (37) 

with the covariantly chiral superfield Z as the Lagrange multiplier superfield. Varying 
the action (37) with respect to Z gives back the original action (35) provided that 

F(K) = nZ(K) - V(Z(K)) (38) 

where the function Z(1Z) is defined by inverting the function 

TZ = V\Z) (39) 

Equations (38) and (39) define the superfield Legendre transform, and imply 

F'{TZ) = Z{TZ) and F"(TZ) = Z'(K) = yj^j^ (40) 

where V" = d 2 V/dZ 2 . The second formula (40) is the duality relation between the 
supergravitational function F and the chiral superpotential V. 

A supersymmetric (local) Weyl transform of the acton (37) can be done entirely in 
superspace. In terms of the field components, the super- Weyl transform amounts to a 
Weyl transform, a chiral rotation and a (superconformal) S'-supersymmetry transfor- 
mation [46]. The chiral density superfield £ appears to be the chiral compensator of 
the super- Weyl transformations, 

£ ->• e m £ (41) 

whose parameter $ is an arbitrary covariantly chiral superfield, V.$ = 0. Under the 
transformation (41) the covariantly chiral superfield 1Z transforms as 

K -)• e- 2<E> (k - |V 2 ) e* (42) 



The super- Weyl chiral superfield parameter $ can be traded for the chiral Lagrange 
multiplier Z by using a generic gauge condition 

Z = Z($) (43) 

where Z(&) is a holomorphic function of $. It results in the action 

= j d*xd 4 6 E^e*** [Z($) + H.c] - J d 4 xd 2 6 £e^V{Z(§)) + H.c. (44) 

Equation (44) has the standard form of the action of a chiral matter superfield 
coupled to supergravity, 



£[$,$]= J d 4 xd A 9E- l n(^,^) 



d 4 xd 2 e£P(^) + R.c. 



(45) 
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in terms of the non-chiral potential fi(3>, $) and the chiral superpotential P($). In our 
case (44) we find 

$) = e $+ * [Z($) + Z($)] , P($) = -e 3 * (46) 
The Kahler potential A'($, <5) is given by 

K = -3\n(--) or tt = -3e- K/3 (47) 
3 

so that the action (45) is invariant under the supersymmetric (local) Kahler- Weyl trans- 
formations 

K($, $) <5) + A($) + A($) , P($) ->■ - e - A(<E>) P($) (48) 

with the chiral superfield parameter A($). It follows that 

£ ->• e A ^£ (49) 

The scalar potential in terms of the usual fields is given by the standard formula 
[47] 

V(0,0) = e*{|!£ + HP| 2 -3|P| 2 }| (50) 

where all the superfields are restricted to their leading field components, $| = (j)(x), 
and we have introduced the notation 



dP dK „ 2 
1 P 



P $ P| 2 = DzP{K-l)D*P (51) 



with = d 2 K/d&d&. Equation (50) can be simplified by making use of the 

Kahler- Weyl invariance (48) that allows one to choose a gauge 

P = 1 (52) 

It is equivalent to the well known fact that the scalar potential (50) is actually governed 
by the single (Kahler- Weyl-invariant) potential 

= fi + ln|P| 2 (53) 

In our case (46) we find 

G = e* + * [Z($) + £($)] + 3($ + <5) + ln(y(Z($)) + ln(V(J(<5)) (54) 
So let us choose a gauge by the condition 

3$ + ln(y(Z($)) = or = e" 3 * (55) 
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that is equivalent to eq. (52). Then the G-potential (54) gets simplified to 



G = e* + * [Z($)+Z(Q)] 



(56) 



There is the correspondence between a holomorphic function F(1Z) in the supergravity 
action (35) and a holomorphic function Z(<&) defining the scalar potential (50), 



To the end of this section, I would like to comment on the standard way of the 
inflationary model building by a choice of $) and -P($) — see eg., refs. [48, 49] 
for a review. 

The factor exp(fT/M| 1 ) in the F-type scalar potential (50) of the chiral matter- 
coupled supergravity, in the case of the canonical Kahler potential, K oc $$, results 
in the scalar potential V oc exp(|$| 2 /Mp^) that is too steep to support chaotic inflation. 
Actually, it also implies r\ ps 1 or, equivalently, M? flaton ps Vq/M^ ps i7 2 . It is known 
as the rj-problem in supergravity [50]. 

As is clear from our discussion above, the ^-problem is not really a supergravity 
problem, but it is the problem associated with the choice of the canonical Kahler po- 
tential for an inflaton superfield. The Kahler potential in supergravity is a (Kahler) 
gauge-dependent quantity, and its quantum renormalization is not under control. Un- 
like the one-field inflationary models, a generic Kahler potential is a function of at least 
two real scalar fields, so it implies a nonvanishing curvature in the target space of the 
non-linear sigma-model associated with the Kahler kinetic term. 3 Hence, a generic 
Kahler potential cannot be brought to the canonical form by a field redefinion. 

To solve the ^-problem associated with the simplest (naive) choice of the Kahler 
potential, on may assume that the Kahler potential K possesses some shift symme- 
tries (leading to its flat directions), and then choose inflaton in one such flat direction 
[52]. However, in order to get inflation that way, one also has to add ("by hand") 
the proper inflaton superpotential breaking the initially introduced shift symmetry, and 
then stabilize the inflationary trajectory with the help of yet another matter superfield. 

The possible alternative is the D-term mechanism [53], where the inflaton particle 
belongs to the matter gauge sector and, as a result, inflation is highly sensitive to gauge 
charges [53]. 

It is worth mentioning that in the (perturbative) superstring cosmology one gets the 
Kahler potential (see e.g., refs. [54, 55]) 



over a Calabi-Yau (CY) space in the type-IIB superstring compactification, thus avoid- 
ing the ^-problem but leading to a plenty of choices ("embarrassment of riches") in the 
String Landscape. 

3 See eg., ref. [51] for more about the non-linear sigma-models. 




(57) 



in the classically equivalent scalar-tensor supergravity. 



K oc log(moduli polynomial)cY 



(58) 
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Finally, one still has to accomplish stability of a given inflationary model in super- 
gravity against quantum corrections. Such corrections can easily spoil the flatness of 
the inflaton potential. The Kahler kinetic term is not protected against quantum correc- 
tions, because it is given by a full superspace integral (unlike the chiral superpotential 
term). The F(1Z) supergravity action (35) is given by a chiral superspace integral, so 
that it is protected against the quantum corrections given by full superspace integrals. 

To conclude this section, we claim that an N = 1 locally supersymmetric extension 
of f(R) gravity is possible. It is non-trivial because the auxiliary freedom has to be 
preserved. The new supergravity action (35) is classically equivalent to the standard 
N — 1 Poincare supergravity coupled to a dynamical chiral matter superfield, whose 
Kahler potential and the superpotential are dictated by a single holomorphic function. 
Inflaton can be identified with the real scalar field component of that chiral matter 
superfield originating from the supervielbein, and thus has the geometrical origin. 

It is worth noticing that the action (35) allows a natural extension in chiral curved 
superspace, due to the last equation (32), namely, 



where W a/ 3 7 is the iV = 1 covariantly-chiral Weyl superfield of the N = 1 superspace 
supergravity, and W 2 = W Q/ 3 7 >V a/37 . The action (59) also has the auxiliary freedom. 
In Supersring Theory, the Weyl-tensor-dependence of the gravitational effective action 
is unambigously determined by the superstring scattering amplitudes or by the super- 
Weyl invariance of the corresponding non-linear sigma-model (see eg., ref. [51]). 

A possible connection of F(1Z) supergravity to the Loop Quantum Gravity was 
investigated in ref. [3]. 



In this section we would like to investigate a possibility of spontaneous supersymmetry 
breaking, without fine tuning, by imposing the condition of the vanishing scalar po- 
tential. Those no-scale supergravities are the starting point of many phenomenological 
applications of supergravity in HEP and inflationary theory, including string theory 
applications — see eg., refs. [56, 57] and references therein. 

The no-scale supergravity arises by demanding the scalar potential (50) to vanish. 
It results in the vanishing cosmological constant without fine-tuning [58]. The no- 
scale supergravity potential G has to obey the non-linear 2nd-order partial differential 
equation, which follows from eq. (57), 




(59) 



7 No-scale F(1Z) Supergravity 



d 2 G dGOG 
<9$<9<S ~ d$ M 
A gravitino mass m 3 / 2 is given by the vacuum expectation value [41] 



(60) 



m 3/2 = (e G ^) 



(61) 
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so that the spontaneous supersymmetry breaking scale can be chosen at will. 
The well known exact solution to eq. (60) is given by 

G = -31og($ + $) (62) 

In the recent literature, the no-scale solution (62) is usually modified by other terms, 
in order to achieve the universe with a positive cosmological constant — see e.g., the 
KKLT mechanism [59]. 

To appreciate the difference between the standard no-scale supergravity solution 
and our 'modified' supergravity, it is worth noticing that the Ansatz (62) is not favoured 
by our potential (56). In our case, demanding eq. (60) gives rise to the lst-order non- 
linear partial differential equation 



3 (e*X' + e $ X') 
where we have introduced the notation 



(63) 



Z(<&) = e -*X(Q) , X ' = % (64) 

in order to get the differential equation in its most symmetric and concise form. 
Accordingly, the gravitino mass (61) is given by 

m 3/2 = ^exp X - [fX + e^x) \ (65) 

I am not aware of any non-trivial holomorphic exact solution to eq. (63). However, 
should it obey a holomorphic differential equation of the form 

X' = e*g{X,$) (66) 



with a holomorphic function g(X, $), eq. (63) gives rise to the functional equation 



3(<7 



e*g + X 



(67) 



Being restricted to the real variables $ = $ = y and X = X = x, eq. (63) reads 

(1 T 

6x> = e y (x' + xf , where x' = — (68) 

dy 

This equation can be integrated after a change of variables, 

x = e~ y u , (69) 
and it leads to a quadratic equation with respect to u' = du/dy, 

(u'f - 6u' + 6u = (70) 

It follows 
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8 Fields from Superfields in F(1Z) Supergravity 

For simplicity, we now set all fermionic fields to zero, and keep only bosonic field 
components of the superfields. It greatly simplies all equations but makes supersym- 
metry to be manifestly broken. Of course, SUSY is restored after adding back all the 
fermionic terms. 

Applying the standard superspace chiral density formula [40, 41, 42] 

J d 4 xd 2 eSC = J d 4 xe{C last + BC Rrst } (72) 

to the action (35) yields its bosonic part in the form 

(-g)- 1/2 L hos = f(R, R; X, X) = F\X) [\R* + AXX] + 3XF(X) + H.c. (73) 

where the primes denote differentiation with respect to the given argument. We have 
used the notation 

X = \B and R* = R + l -e abcd R abcd = R + iR (74) 

The R does not vanish in F(1Z) supergravity, and it represents the axion field that is 
the pseudo-scalar superpartner of real scalaron field in our construction. 
Varying eq. (73) with respect to the auxiliary fields X and X, 

dL^os _ <9Lbos _ Q ,rjc\ 

~dx~~~dT~ { } 

gives rise to the algebraic equations on the auxiliary fields, 

3F + X(AF' + IF 1 ) + AXXF" + ±F"R* = (76) 

and its conjugate 

3F + X(AF' + IF') + AXXF" + $F"R* = (77) 

where F = F(X) and F = F(X). The algebraic equations (76) and (77) cannot be 
explicitly solved for X in a generic F(1Z) supergravity. 

To recover the standard (pure) supergravity in our approach, let us consider the 
simple special case when 

F" = or, equivalents, F(K) = f - \f{R, (78) 

with some complex constants f and fi, where Re fx > 0. Then eq. (76) is easily 
solved as 

5(Re/i) 
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Substituting this solution back into the Lagrangian (73) yields 

L = -\{Reh)R + = -\Ml x R - A (80) 

where we have introduced the reduced Planck mass M Pi , and the cosmological con- 
stant A as 

Re/i = §M|! and A = (8 i) 
It is the standard pure supergravity with a negative cosmological constant [40, 41, 42]. 



9 Generic 7Z 2 supergravity, and AdS Bound 

The simplest non-trivial F(1Z) supergravity is obtained by choosing F" = const. ^ 
that leads to the 1Z 2 - supergravity defined by a generic quadratic polynomial in terms 
of the scalar supercurvature [8]. 

Let us recall that the stability conditions in f(R) -gravity are given by eqs. (12) in 
the notation (9). In the notation (73) used here, ie. when f(R) = —\Mp l f(R), one 
gets the opposite signs, 

f'(R)<0 (82) 

and 

f"(R) > (83) 

The first (classical stability) condition (82) is related to the sign factor in front of 
the Einstein-Hilbert term (linear in R) in the f(R) -gravity action, and it ensures that 
graviton is not a ghost. The second (quantum stability) condition (83) guarantees that 
scalaron is not a tachyon. 

Being mainly interested in the inflaton part of the bosonic f(R) -gravity action that 
follows from eq. (73), we set both gravitino and axion to zero, which also implies 
i?* = R and a real X. 

In F(R) supergravity the stability condition (82) is to be replaced by a stronger 
condition, 

F'(X) < (84) 

It is easy to verify that eq. (82) follows from eq. (84) because of eq. (75). Equation 
(84) also ensures the classical stability of the bosonic f(R) gravity embedding into the 
full F{1Z) supergravity against small fluctuations of the axion field. 

In this Section we investigate the most general Ansatz (with F" = const. ^ 0) 
that leads to the simplest non-trivial toy-model of F(R) supergravity with the master 
function 

F(n) = h- l -hn+ l -f 2 n 2 (85) 
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having three coupling constants / , fi and f 2 . We take all of them to be real, since 
we ignore this potential source of CP-violation here (see, however, the Outlook in 
Sec. 15). As regards the mass dimensions of the quantities introduced, we have 

[F] = [f ] = 3, [#] = [/!] =2, and [R] = [f 2 ] = 1 (86) 

The bosonic Lagrangian (73) with the function (85) reads 

(-g)- 1/2 L hos = llf 2 X 3 - 7frX 2 + (lf 2 R + 6/ ) X - \f\R (87) 

Hence, the auxiliary field equation (75) takes the form of a quadratic equation, 

f f 2 X 2 - lj x X + \Rj 2 + 3/ = (88) 

whose solution is given by 

7 



X. 



± 



3-11 



(89) 



where we have introduced the maximal scalar curvature 

*-= 2^1 - 3 1 (90) 

Equation (89) obviously implies the automatic bound on the scalar curvature (from one 
side only). In our notation, it corresponds to the (AdS) bound on the scalar curvature 
from above, 

R < Pmax (91) 

The existence of the built-in maximal (upper) scalar curvature (or the AdS bound) 
is a nice bonus of our construction. It is similar to the factor a/1 — v 2 /c 2 in Special 
Relativity. Yet another close analogy comes from the Born-Infeld non-linear extension 
of Maxwell electrodynamics, whose (dual) Hamiltonian is proportional to [51] 

(i - \l i - ^/^Lx - ^7#Lx + (e x W/ELJiL ) (92) 

— * — * 

in terms of the electric and magnetic fields E and H, respectively, with their maximal 
values. For instance, in String Theory one has E mSLX = H mSLX = (27ia')~ 1 [51]. 

Substituting the solution (89) back into eq. (87) yields the corresponding f(R)- 
gravity Lagrangian 

f(mJ-lMl 2 " 73 8. 

J±{ ] 11 h 3 3 • ll 2 /I 

(93) 



19 h R T\ Hr(%J2){R^-R) 



32 . n ^ 1 y 11 \ 3 3 
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Expanding eq. (93) into power series of R yields 

f ± (R) = -A ± - a± R + b±R 2 + 0(R 3 ) (94) 
whose coefficients are given by 

A± = 3^-/1 (i?max " J-hn^) ± (| /2 ) ^ (95) 



and 



3 2 -11 V 11 V 32 • 

h 



19 / 2 / 2 

Oi = 7T A =F \ / TT^max 7^/2 ) (96) 



(97) 



>± 



1 li? max y2 • 3 2 

Those equations greatly simplify when f = 0. One finds [5, 8] 
where we have chosen 



(98) 



/1 = ^pi (99) 

in order to get the standard normalization of the Einstein-Hilbert term that is linear 
in R. Then, in the limit R m3X — > +00, both functions f±\R) reproduce General 
Relativity. In another limit R -» 0, one finds a vanishing or positive cosmological 
constant, 

A W = and Af = A-^M^ max (100) 

The stability conditions are given by eqs. (82), (83) and (84), while the 3rd condi- 
tion implies the 2nd one. In our case (93) we have 

f±(R) = - 327^/1 ± (f /2 ) ^^max - R < (101) 

and 



u R) = T d) > (102) 



while eqs. (84), (85) and (89) yield 



±y^(iw-#) <7^| ( 103 ) 

It follows from eq. (102) that 

/ 2 (+) <0 and ft ] >0 (104) 
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Then the stability condition (83) is obeyed for any value of R. 

As regards the (-)-case, there are two possibilities depending upon the sign of fx. 
Should fi be positive, all the remaining stability conditions are automatically satisfied, 
ie. in the case of both f ( 2 ) > Oand/J^ > 0. 

Should fi be negative, f[~^ < 0, we find that the remaining stability conditions 
(101) and (103) are the same, as they should, while they are both given by 

r r 19 ft _ _ 3-5 fl 2 / ins nn^ 

K < K ^ 2 3 . n /a 2 - 2 3 . H J2 /2 " ^max U^j 

As regards the (+)-case, eq. (103) implies that fi should be negative, f\ < 0, 
whereas then eqs. (101) and (103) result in the same condition (105) again. 

Since R™l x < -Rmax, our results imply that the instability happens before R reaches 
-Rmax in all cases with negative f\ . 

As regards the particularly simple case (98), the stability conditions allow us to 
choose the lower sign only. 

A different example arises with a negative f\. When choosing the lower sign (ie. a 
positive f 2 ) for definiteness, we find 



f-(R) = -^fo 



h 



2-7 3 
3 3 • ll 2 



ft 



ft 



+ 3^I 1/11 R+ STi (P 2 ) ^~ R f' 2 



(106) 



Demanding the standard normalization of the Einstein-Hilbert term in this case implies 

_3 4 -ll/M| 1 19 V 

i? max - -^jr \~y + ^tyy l/il ) (107) 

where we have used eq. (96). It is easy to verify by using eq. (95) that the cosmological 
constant is always negative in this case, and the instability bound (105) is given by 

3 4 -llMg 1 /Mj, 19_LA 
2 3 fl V 22 3 2 • 11 



Rmlx - ( + o9 ~i ) < ^max (108) 



The f-(R) function of eq. (93) can be rewritten to the form 



f( R ) ~ 33 . H2^~ 3 2 . n A i? max- 32 ; ^f\R + h\j ^-^j{R max - Rf /2 (109) 

where we have used eq. (90). There are three physically different regimes: 

(i) the high-curvature regime, R < and \R\ 3> R max . Then eq. (109) implies 

f(R) » -A h -a h R + c h \R\ 3/2 (110) 
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whose coefficients are given by 



a _ 2-7 7 3 A 3 

A h — 09 1 1 h K 



3 2 -ll J1 max 3 3 • ll 2 /| ' 

19 
3 2 • 11 1 



!9 , 

a h = ^r—fi , (ill) 



Cft - Vn (l /2 

(ii) the low-curvature regime, \R/R max \ <C 1. Then eq. (109) implies 

/(i?) » -A, - ai R , (112) 

whose coefficients are given by 



(113) 

where we have used eq. (96). 

(iii) the near-the-bound regime (assuming that no instability happens before it), 
R = ^max + SR, 5R < 0, and \5R/R mBX \ < 1. Then eq. (109) implies 

f(R) « -A 6 + a fe + c 6 \5R\ 3/2 (1 14) 

whose coefficients are 

1 7 3 f 3 

iv £> — g Jl-fl-max 33 . J| ' 

a 6 = a h , (115) 
Cb = {fl (l /2 

The cosmological dynamics may be either directly derived from the gravitational 
equations of motion in the /(i?) -gravity with a given function f(R), or just read off 
from the form of the corresponding scalar potential of a scalaron (see below). For 
instance, as was demonstrated in ref. [5] for the special case / = 0, a cosmological 
expansion is possible in the regime (i) towards the regime (ii), and then, perhaps, to 
the regime (iii) unless an instability occurs. 

However, one should be careful since our toy-model (85) does not pretend to be 
viable in the low-curvature regime, eg., for the present Universe. Nevertheless, if one 
wants to give it some physical meaning there, by identifying it with General Relativity, 
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then one should also fine-tune the cosmological constant in eq. (1 13) to be "small" 
and positive. We find that it amounts to 



3 4 -7 2 -llM^ 



pi 



25-192 /I 



R 



A=0 



(116) 



with the actual value of i? max to be "slightly" above of that bound, R 



max > R \ „• 



It 



is also posssible to have the vanishing cosmological constant, A; = 0, when choosing 
= -R A=n . It is worth mentioning that it relates the values of -R max and f 2 . 



Rm&y 

The particular 7?. 2 -supergravity model (with f = 0) was introduced in ref. [5] in 
an attempt to get viable embedding of the Starobinsky model into F(7£)-supergravity. 
However, it failed because, as was found in ref. [5], the higher-order curvature terms 
cannot be ignored in eq. (98), ie. the i? n -terms with n > 3 are not small enough 
against the i? 2 -term. In fact, the possibility of destabilizing the Starobinsky inflationary 
scenario by the terms with higher powers of the scalar curvature, in the context of f(R) 
gravity, was noticed earlier in refs. [60, 61]. The most general Ansatz (85), which is 
merely quadratic in the supercurvature, does not help for that purpose either. 

For example, the full /(i?) -gravity function f-(R) in eq. (98), which we derived 
from our 7?. 2 -supergravity, gives rise to the inflaton scalar potential 



V(y) = V (Ue v + 3) {e~ y - l) 2 
where V = (3 3 /2 6 )MpJ /f . The corresponding inflationary parameters 



(117) 



and 



e(y) = 3 



v(y) = 



[liey + 3)(e-v -1) 



1 

> - 

- 3 



2 (lief + he- y + I2e~ 2y ) 2 

3 (lief + 3)(e-f-l) 2 ~ 3 



(118) 



(119) 



are not small enough for matching the WMAP observational data. A solution to this 
problem is given in the next Sec. 10. 



10 Chaotic inflation in F(1Z) Supergravity 

Let us take now one more step further and consider a new Ansatz for F(1Z) function 
in the cubic form 

F(K) = -hfiK + \hK 2 - IW 3 (120) 

whose real (positive) coupling constants /i,2,3 are of (mass) dimension 2, 1 and 0, 
respectively. Our conditions on the coefficients are 

h » i , fi » h (i2i) 
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The first condition is needed to have inflation at the curvatures much less than Mp } 
(and to meet observations), while the second condition is needed to have the scalaron 
(inflaton) mass be much less than Mpi, in order to avoid large (gravitational) quantum 
loop corrections after the end of inflation up to the present time. 

The bosonic action is given by eq. (73). In the case of a real scalaron it reduces to 



Lj^Tg = IF 1 



-R + AX 2 
3 



+ 6XF 



(122) 



so that the real auxiliary field is a solution to the algebraic equation 

'1 



3F + llF'X + F" 



-R + 4X Z 







(123) 



Stability of the bosonic embedding in supergravity requires F'(X) < (Sec. 9). 
In the case (120) it gives rise to the condition /| < fif 3 . For simplicity here, we will 
assume a stronger condition, 

ft « hh (124) 

Then the second term on the right-hand- side of eq. (120) will not affect inflation, as is 
shown below. 

Equation (122) with the Ansatz (120) reads 



L = -5f 3 X 4 + llf 2 X 3 - (7/x + \hR)X 2 + lf 2 RX - \f,R 
and gives rise to a cubic equation on X, 



\20fJ \10f 3 30 J 



30/ ; 



R = 



(125) 



(126) 



We find three consecutive (overlapping) regimes. 

• The high curvature regime including inflation is given by 



SR < and 



Rn 



> 



fi 



hh 



1/3 



(127) 



where we have introduced the notation R = 21/i// 3 > and 5R = R + 
R . With our sign conventions we have R < during the de Sitter and matter 
dominated stages. In the regime (127) the / 2 -dependent terms in eqs. (125) and 
(126) can be neglected, and we get 



X 2 - -^5R 



and 



L=J i R+ ^ R+R ° 



(128) 
(129) 
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It closely reproduces the Starobinsly inflationary model (Sec. 2) since inflation 
occurs at \R\ 3> Rq. In particular, we can identify 



f - 15M " mm 

It is worth mentioning that we cannot simply set f 2 = in eq. (120) because it 
would imply X = and L = —^-R for 5R > 0. As a result of that the scalar 
degree of freedom would disappear that would lead to the breaking of a regular 
Cauchy evolution. Therefore, the second term in eq. (120) is needed to remove 
that degeneracy. 

• The intermediate (post-inflationary) regime is given by 

\5R 
Ro 

In this case X is given by a root of the cubic equation 



<1 (131) 



30X 3 + (6R)X + = (132) 

J3 

It also implies that the 2nd term in eq. (126) is always small. Equation (132) 
reduces to eq. (128) under the conditions (127). 

The low-curvature regime (up to R = 0) is given by 

SR f f 2 \ ^ 3 
SR > o and — > -p- (133) 

Ro V/l/3/ 

It yields 

x ^ mhr) (134) 

and 

It is now clear that fi should be equal to 3Mp t /2 in order to obtain the correctly 
normalized Einstein gravity at \R\ R . In this regime the scalaron mass 
squared is given by 

1 fsRoM^ 21/, „ 63M|, 



3|/"(i2)| 4/| " 4/| F1 8/. 



2 



in agreement with the case of the absence of the 1Z 3 term, studied in the previous 
section. The scalaron mass squared (136) is much less than Mp t indeed, due to 
the second inequality in eq. (121), but it is much more than one at the end of 
inflation (~ M 2 ). 



28 



It is worth noticing that the corrections to the Einstein action in eqs. (129) and (135) 
are of ther same order (and small) at the borders of the intermediate region (131). 

The roots of the cubic equation (126) are given by the textbook (Cardano) formula 
[62], though that formula is not very illuminating in a generic case. The Cardano for- 
mula greatly simplifies in the most interesting (high curvature) regime where inflation 
takes place, and the Cardano discriminant is 

3 

D 



(so) <° ( 137 ) 



It implies that all three roots are real and unequal. The Cardano formula yields the 
roots 

X,, 3 « f cos ( 4f3 / WR/fi + cjj + M (138 ) 

where the constant 61,2,3 takes the values (71-/6, 5n/Q, Stt/2). 

As regards the leading terms, eqs. (125) and (138) result in the (— i?) 3//2 correction 
to the (R + _R 2 )-terms in the effective Lagrangian in the high-curvature regime \R\ ^> 
fi/fi- l n order to verify that this correction does not change our results under the 
conditions (127), let us consider the /(i?) -gravity model with 

f(R) =R- b(-Rf 12 - aR 2 (139) 

whose parameters a > and b > are subject to the conditions a> 1 and b/a 2 <C I. 
It is easy to check that f'(R) > for R e (— 00, 0], as is needed for (classical) stability. 

Any /(it!) gravity model is classically equivalent to the scalar-tensor gravity with 
certain scalar potential (Sec. 3). The scalar potential can be calculated from a given 
function f(R) along the standard lines (Sec. 3). We find (in the high curvature regime) 

K M = _L (l _ e -^ + _^e-*(e»-l) 3 / 2 (.40) 

in terms of the inflaton field y. The first term of this equation is the scalar potential 
associated with the pure (R + R 2 ) model, and the 2nd term is the correction due to the 
i? 3 / 2 -term in eq. (139). It is now clear that for large positive y the vacuum energy in 
the first term dominates and drives inflation until the vacuum energy is compensated 
by the y-dependent terms near e y = 1 . 

It can be verified along the lines of ref. [35] that the formula for scalar perturbations 
remains the same as that for the model (7), ie. A^ ps N 2 M 2 nf / (247r 2 M| 1 ), where iV is 
the number of e-folds from the end of inflation. So, to fit the observational data, one 
has to choose 

h ~ 5iV 2 /(87r 2 A^) ^ 6.5 • 10 10 (iV e /50) 2 (141) 

Here the value of is taken from ref. [12] and the subscript 1Z has a different mean- 
ing from the rest of this review. 
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We conclude that the model (120) with a sufficiently small f 2 obeying the condi- 
tions (121) and (124) gives a viable realization of the chaotic (R + i? 2 )-type inflation 
in supergravity. The only significant difference with respect to the original (R+ R 2 ) 
inflationary model is the scalaron mass that becomes much larger than M in supergrav- 
ity, soon after the end of inflation when 5R becomes positive. However, it only makes 
the scalaron decay faster and creation of the usual matter (reheating) more effective. 

The whole series in powers of 1Z may also be considered, instead of the limited 
Ansatz (120). The only necessary condition for embedding inflation is that / 3 should 
be anomalously large. When the curvature grows, the 7?. 3 -term should become impor- 
tant much earlier than the convergence radius of the whole series without that term. 
Of course, it means that viable inflation may not occur for any function F(1Z) but only 
inside a small region of non-zero measure in the space of all those functions. However, 
the same is true for all known inflationary models, so the very existence of inflation 
has to be taken from the observational data, not from a pure thought. 

The results of this Section can be considered as the viable alternative to the ear- 
lier fundamental proposals [52, 53] for realization of chaotic inflation in supergravity. 
But inflation is not the only target of our construction. As is well known [21, 22, 63], 
the scalaron decays into pairs of particles and anti-particles of quantum matter fields, 
while its decay into gravitons is strongly suppressed [64]. It thus represents the uni- 
versal mechanism of viable reheating after inflation and provides a transition to the 
subsequent hot radiation-dominated stage of the Universe evolution. In its turn, it 
leads to the standard primordial nucleosynthesis (BBN) after. In F(R) supergravity 
the scalaron has a pseudo-scalar superpartner (axion) that may be the source of a strong 
CP-violation and then, subsequently, lepto- and baryo-genesis that naturally lead to 
baryon (matter-antimatter) asymmetry [68, 69] — see Sees. 13 and 15 for more. 

11 Cosmological Constant in F(7Z) Supergravity 

The Standard (A-CDM) Model in cosmology gives a phenomenological description 
of the observed Dark Energy (DE) and Dark Matter (DM). It is based on the use 
of a small positive cosmological constant A and a Cold Dark Matter (CDM), and is 
consistent with all observations coming from the existing cosmological, Solar system 
and ground-based laboratory data. However, the A-CDM Model cannot be the ultimate 
answer to DE, since it implies its time-independence. For example, the 'primordial' 
DE responsible for inflation in the early Universe was different from A and unstable. 
The dynamical (ie. time-dependent) models of DE can be easily constructed by using 
the f(R) gravity theories, defined via replacing the scalar curvature R by a function 
f(R) in the gravitational action. The f(R) gravity provides the self-consistent non- 
trivial alternative to the A-CDM Model. Viable f(R) -gravity -based models of the 
current DE are also known [70, 71, 72], and the combined inflationary-DE models are 
possible too [27]. 

The natural question arises, whether F(1Z) supergravity is also capable to describe 
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the present DE and eg., a positive cosmological constant. It is non-trivial because the 
standard (pure) supergravity can only have a zero or negative cosmological constant. In 
this Section we further extend the Ansatz used in Sec. 10 for the F-function, and apply 
it to get a positive cosmological constant in the regime of low spacetime curvature. 

Throughout this Section we use the units c — h — M P \ = 1. We recall that 
an AdS-spacetime has a positive scalar curvature, and a dS-spacetime has a negative 
scalar curvature in our notation. 

The embedding of f(R) gravity into F(1Z) supergravity is given by (Sec. 8) 

f(R) = f(R,X(R)) (142) 
where the function f(R, X) (or the gravity Lagrangian C) is defined by 

1 



C = f(R,X)=2F'(X) 



-R + AX 2 
3 



+ 6XF(X) (143) 



and the function X = X(R) is determined by solving an algebraic equation, 

T = (144) 

The cosmological constant in F(1Z) supergravity is thus given by 

A=-f(0,X o ) (145) 

where X = X(0). It should be mentioned that X represents the vacuum expectation 
value of the auxiliary field X that determines the scale of the supersymmetry breaking. 
Both inflation and DE imply X ^ 0. 

To describe DE in the present Universe, ie. in the regime of low spacetime cur- 
vature R, the function f(R) should be close to the Einstein-Hilbert (linear) function 
/eh(-R) with a small positive A, 

\f(R) - /EH(i2)| « |/ EH (i2)| , \f'(R) - /ehI « 1, \Rf"(R)\ « 1 (146) 

ie. f(R) —\R — A for small R with the very small and positive A ss 10 _118 (Mp 1 ). 
Equations (143) and (145) imply 

A = -8F'(X )X 2 - 6X F(X ) (147) 

where X is a solution to the algebraic equation 

4X 2 F"(X ) + 11X F'(X ) + 3F(X ) = (148) 

As is clear from eq. (147), to have A ^ 0, one must have X ^ 0, ie. a (spontaneous) 
supersymmetry breaking. However, in order to proceed further, we need a reasonable 
Ansatz for the F-function. 
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The simplest opportunity is given by expanding the function F(TZ) in Taylor series 
with respect to 1Z. Since the N = 1 chiral superfield 1Z has X as its leading field 
component (in ^-expansion), one may expect that the Taylor expansion is a good ap- 
proximation as long as \Xq\ <C l(Mpi). As was demonstrated in Sec. 10, a viable 
(successful) description of inflation is possible in F(1Z) supergravity, when keeping 
the cubic term TZ? in the Taylor expansion of the F(1Z) function. It is, therefore, natu- 
ral to expand the function F up to the cubic term with respect to 1Z, and use it as our 
Ansatz here, 

F{n) = f - \hU + ^f 2 U 2 - l -hK* (149) 

with some real coeffieints f , f 1: f 2) fs- The Ansatz (149) differs from the one used in 
eq. (120) by the presence of the new parameter f only. It is worth emphasizing here 
that f is not a cosmological constant because one still has to eliminate the auxiliary 
field X. The stability conditions (Sec. 9) imply 

A > , f 2 > , h > (150) 

and 

fl < hh (15D 

Inflation requires f 3 3> 1 and /f 3> f\. 4 As was already found in Sec. 10, in order 
to meet the WMAP observations, the parameter / 3 should be approximately equal to 
6.5 • 10 10 (iV e /50) 2 . The cosmological constant in the high-curvature regime does not 
play a significant role and may be ignored there. 

In the low curvature regime, in order to recover the Einstein-Hilbert term, one has 
to fix fi = 3/2 (Sec. 10). Then the Ansatz (149) leads to the gravitational Lagrangian 

f(R, X) = -5/ 3 X 4 + llf 2 X 3 - l -h (r + |£) X 2 + (Vo + \hR^ X- l -R 

(152) 

and the auxiliary field equation 

x3 -w^ + U R+ w) x -w^ R+9f '^ (153 » 

whose formal solution is available via the standard Cardano (Viete) formulae [62]. 
In the low-curvature regime we find a cubic equation for X in the form 

'Linearizing' eq. (154) with respect to X brings the solution X = 2f /7 whose 
substitution into the action (152) gives rise to a negative cosmological constant, A = 
—6/o/7. This way we recover the standard supergravity case. 

4 The stronger condition /| <C /1/3 has been used in Sec. 10 for simplicity. 
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Equations (152) and (154) allow us to write down the exact eq. (145) for the cos- 
mological constant in the factorized form 



A(X 



11/2 



X (X -X_)(X -X + ) 



(155) 



where X± are the roots of the quadratic equation x 2 — jffex + 

21 



0, ie. 



X 4 



22/ 2 



2 3 • 11 



(156) 



Since f f 2 is supposed to be very small, both roots X± are real and positive. 

Equation (155) implies that A > when either (I) X < 0, or (II) X is inside the 
interval (X_,X + ). 

By using Matematica we were able to numerically confirm the existence of solu- 
tions to eq. (154) in the region (I) when f < 0, but not in the region (II). So, to this 
end, we continue with the region (I) only. All real roots of eq. (154) are given by 



(X ) 2 



2 a/— Q cos 
2 a/— Q cos 
2 a/— Q cos 



$\ 11/2 

3/ 20/ 3 



$ + 2tt 
3 

■0 + 47T 



+ 



+ 



II/2 
2O/3 

II/2 
2O/3 



(157) 



in terms of the Cardano-Viete parameters 



Q = 

R = 



II/2 
2 2 • 5/3 2 4 

25.5VI 



7 2 



• 5 2 /| 
3/o 
2 2 • 5/3 



and the angle 1? defined by 



II/2 
2O/3 ' 

H 3 /l 
2 6 • 53/1 



cos-# = 



1 



20/3 v"Y g + 3/o 



(158) 



R 



(159) 



The Cardano discriminant reads D = R 2 + Q 3 . All three roots are real provided that 
D < 0. It is known to be the case in the high-curvature regime (Sec. 10), and it is also 
the case when / is extremely small. Under our requirements on the parameters the 
angle d is very close to zero, so the relevant solutions X < are given by the 2nd and 
3rd lines of eq. (157), with X ps / /10. 

We thus showed that it is possible to have a positive cosmological constant (at 
low spacetime curvature) in the particular F(1Z) supergravity (without its coupling to 
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super-matter) described by the Ansatz (149). The same Ansatz is applicable for de- 
scribing a viable chaotic inflation in supergravity (at high spacetime curvature). The 
positive cosmological constant was achieved as the non-linear effect (with respect to 
the superspace curvature and spacetime curvature) in the relatively narrow part of the 
parameter space (it is, therefore, highly constrained). It also implies the apparent vio- 
lation of the Strong Energy Condition in our model. 

Of course, describing the DE in the present Universe requires enormous fine-tuning 
of our parameters in the F-function. However, it is the common feature of all known 
approaches to the DE. Our analysis does not contribute to 'explaining' the smallness 
of the cosmological constant. 

12 Nonminimal Scalar- Curvature Coupling in Gravity 
and Supergravity, and Higgs inflation 

One may pursue different strategies in a theoretical search for inflaton. For instance, 
inflaton may be either a new exotic particle or something that we already know 'just 
around the corner'. In this review I advocate the second ('economical') approach. 
Perhaps, the most 'economical' possibilities are given by Higgs inflation [73, 74, 76] 
and Starobinsky (R + R 2 ) inflation (Sec. 2). 

According to the cosmology textbooks, a Higgs particle of the Standard Model 
cannot serve as inflaton because the SM parameters are A p=* 1, m H ~ 10 2 GeV, and 
(ST/T) 1, whereas inflation requires (see Sec. 4) A ps 10~ 13 , m inf ps 10 13 GeV, 
and (ST/T) ps 1CT 5 . Nevertheless, it is possible to reach the required values when 
assuming that Higgs particle is nonminimally coupled to gravity [73, 74, 76]. For 
instance, adding the nonminimal coupling of the Higgs field to the scalar spacetime 
curvature is natural in curved spacetime because it is required by renormalization [77]. 
It is worth recalling here that in the Starobinsky inflaton [21, 22] the inflaton field is 
the spin-0 part of spacetime metric. 

In this Section we compare the inflationary scalar potential, derived by the use of 
the nonminimal coupling [73, 74, 76], with the scalar potential that follows from the 
(R + R 2 ) inflationary model (Sec. 2), and confirm that they are the same. Then we 
also upgrade that equivalence to supergravity. In this section we set M P1 = 1 too. 

The original motivation of refs. [73, 74, 76] is based on the assumption that there 
is no new physics beyond the Standard Model up to the Planck scale. Then it is natural 
to search for the most economical mechanism of inflation by identifying inflaton with 
Higgs particle. Later, in this Section, we assume that there is the new physics beyond 
the Standard Model, and it is given by s upersymmetry. Then it is quite natural to 
search for the most economical mechanism of inflation in the context of supergravity. 
Moreover, we do not have to identify our inflaton with a Higgs particle of the Minimal 
Super symmetric Standard Model. 
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Let us begin with the 4D Lagrangian 

£j = V=9~j + e0j)^J + - ^(0j)} (160) 

where we have introduced the real scalar field (f>j(x), nonminimally coupled to gravity 
(with the coupling constant £) in Jordan frame, with the Higgs-like scalar potential 

V{h) = \{ti-v 2 ) 2 (161) 

The action (160) can be rewritten to Einstein frame by redefining the metric via a 
Weyl transformation, 

< r = jrfm <162) 

It gives rise to the standard Einstein-Hilbert term (—\R) for gravity in the Lagrangian. 
However, it also leads to a nonminimal (or noncanonical) kinetic term of the scalar field 
0j. To get the canonical kinetic term, a scalar field redefinition is needed, 0j — > y?(</>j), 
subject to the condition 

(163) 

As a result, the non-minimal theory (160) is classically equivalent to the standard 
(canonical) theory of the scalar field ip(x) minimally coupled to gravity, 

£e = V=9 {-\r + \g^d^d v f - 1%) J (164) 
with the scalar potential 

[l+^j(W] 2 

Given a large positive £ >> 1, in the small field limit one finds from eq. (163) that 
pa if, whereas in the large f 1 imit one gets 

¥>«y|log(l + £$) (166) 



#J 1 + 



Then eq. (165) yields the scalar potential: 
(i) in the very small field limit, y? < ^/§£~\ as 



Ks(y?) ~ \^ (167) 



(ii) in the small field limit, * /|£ 1 < f <C \ /§, as 



2' 

A 



~ 7^V% (168) 



6£ 
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(iii) and in the large field limit, ip 3> y §£ 1 , as 



3* 



(169) 



We have assumed here that £ 3> 1 and t)^ C 1. 

Identifying inflaton with Higgs particle requires the parameter v to be of the order 
of weak scale, and the coupling A to be the Higgs boson selfcoupling at the inflationary 
scale. The scalar potential (169) is perfectly suitable to support a slow -roll inflation, 
while its consistency with the WMAP normalization condition (Sec. 4) for the ob- 
served CMB amplitude of density perturbations at the e-foldings number iV e = 55 
gives rise to the relation £/v^A w 5 ■ 10 4 [73, 74, 76]. 

The scalar potential (168) corresponds to the post- inflationary matter-dominated 
epoch described by the oscillating inflaton field ip with the frequency 

u> = yir 1 = M inf (170) 

When gravity is extended to 4D, N = 1 supergravity, any physical real scalar field 
should be complexified, becoming the leading complex scalar field component of a 
chiral (scalar) matter supermultiplet. In a curved superspace of N = 1 supergravity, the 
chiral matter supermultiplet is described by a covariantly chiral superfield $ obeying 
the constaraint = 0. The standard (generic and minimally coupled) matter- 

supergravity action is given by in superspace by eqs. (45) and and (47), namely, 



'MSG 



= -zj d A xd A 6E~ 1 exp ~K($,$) d 4 xd 2 0EW($) + H.c.j (171) 

in terms of the Kahler potential K = — 3 log(— and the superpotential W of the 
chiral supermatter, and the full density E and the chiral density S of the superspace 
supergravity (Sec. 5). 

The non-minimal matter- supergravity coupling in superspace reads 

S NM = J d A xd 2 d£X(<&)K + H.c. (172) 

in terms of the chiral function and the N=l chiral scalar supercurvature super- 

field 1Z obeying V»7£ = 0. In terms of the field components of the superfields the 
non-minimal action (172) is given by 

J d 4 xd 2 9£X($)n + H.c. = ~J d A xyf=gX{4> c )R + H.c. + . . . (173) 

stand for the fermionic terms, and C = $| = 4> + ix is the leading complex scalar field 
component of the superfield $. Given X($) = — £$ 2 with the real coupling constant 
£, we find the bosonic contribution 

•Wos. = \i I d A x^T~g (0 2 - x 2 ) R (174) 
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It is worth noticing that the supersymmetrizable (bosonic) non-minimal coupling reads 
# + (^) 2 ] R, not {^ C )R. 

Let us now introduce the manifestly supersymmetric nonminimal action (in Jordan 
frame) as 

S = Smsg + •S'nm (175) 
In curved superspace of N = 1 supergravity the (Siegel's) chiral integration rule 



d 4 xd 2 6£C ch 



d A xd A 6E 



ch 



K 



(176) 



applies to any chiral superfield Lagrangian £ ch with V«£ c h = 0. It is, therefore, 
possible to rewrite eq. (172) to the equivalent form 



S NM = J d A xd A 9E' 1 [X($)+X($j\ 



(177) 



We conclude that adding S'nm to Smsg is equivalent to the simple change of the Q- 
potential as (cf. ref. [78]) 



n q nm = n + x($) + x($) 

It amounts to the change of the Kahler potential as 



(178) 



K 



NM 



-3 In 



_ K/S _ X($)+X($) 
3 



(179) 



The scalar potential in the matter-coupled supergravity (171) is given by eq. (57), 



v{<t>A) 



d 2 G \ 1 dGdG ' 



d<pd(f) J dcj) d(j> 
in terms of the Kahler-gauge-invariant function (53), ie 

G = K + \n\W\ 2 
Hence, in the nonminimal case (175) we have 



(180) 



G 



NM 



i^NM +ln|W| 



(181) 



(182) 



Contrary to the bosonic case, one gets a nontrivial Kahler potential Xnm, ie- a 
Non-Linear Sigma-Model (NLSM) as the kinetic term of C = + i\ (see ref. [51] for 
more about the NLSM). Since the NLSM target space in general has a nonvanishing 
curvature, no field redefinition generically exist that could bring the kinetic term to the 
free (canonical) form with its Kahler potential Kf ree = $$. 
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Let's now consider the full action (175) under the slow-roll condition, ie. when the 
contribution of the kinetic term is negligible. Then eq. (175) takes the truly chiral form 

S ch . = J d 4 xd 2 9S [X($)K + W($)] + H.c. (183) 

When choosing X as the independent chiral superfield, S'ch. can be rewritten to the 
form 

S ch . = J d A xd 2 9£ [XII - Z{X)\ + H.c. (184) 

where we have introduced the notation 

Z{X) = -W($(X)) (185) 

In its turn, the action (184) is equivalent to the chiral F(1Z) supergravity action 
(35), whose function F is related to the function Z via Legendre transformation (Sec. 6) 

Z = X1Z — F , F\K) = X and Z\X) = TZ (186) 

It implies the equivalence between the reduced action (183) and the corresponding 
F(1Z) supergravity whose F-finction obeys eq. (186). 

Next, let us consider the special case of eq. (183) when the superpotential is given 

by 

= ^m$ 2 + iA$ 3 (187) 

with the real coupling constants m > and A > 0. The model (187) is known as the 
Wess-Zumino (WZ) model in 4D, N — 1 rigid supersymmetry. It has the most general 
renormalizable scalar superpotential in the absence of supergravity. In terms of the 
field components, it gives rise to the Higgs-like scalar potential. 
For simplicity, let us take a cubic superpotential, 

W 3 (<S>) = ^A$ 3 (188) 

or just assume that this term dominates in the superpotential (187), and choose the 
X($)-function in eq. (183) in the form 

X($) = -£$ 2 (189) 

with a large positive coefficient £, £ > and £ 3> 1, in accordance with eqs. (173) and 
(174). 

Let us also simplify the F-function of eq. (120) by keeping only the most relevant 
cubic term, 

F 3 (n) = -IfzK 3 (190) 
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It is straightforward to calculate the Z-function for the F-function (190) by using 
eq. (186). We find 



-X = -f 3 7Z 2 and Z'(X) = J-— (191) 
Integrating the last equation with respect to X yields 

' ^ 3 J 3 

where we have used eq. (189). In accordance to eq. (185), the F(7£.)- supergravity 
Z -potential (192) implies the superpotential 

2\/2 f 3 / 2 

J 3 

It coincides with the superpotential (188) of the WZ-model, provided that we identify 
the couplings as 

/a = S£ (194) 

A z 

We conclude that the original nonminimally coupled matter- supergravity theory 
(175) in the slow-roll approximation with the superpotential (188) is classically equiv- 
alent to the F (71) -supergravity theory with the F-function given by eq. (190) when 
the couplings are related by eq. (194). 

The inflaton mass M in the supersymmertic case, according to eqs. (130) and (194), 
is given by 

*£, = g 

Since the value of M in f is fixed by the WMAP normalization (Sec. 4), the value of £ in 
the supersymmetric case is £g = (45/32)^ os , or £ susy 10 3 , ie. is lower than that 
in the bosonic case. We have asssumed here that A w 0(1). 

The established equivalence begs for a fundamental reason. In the high-curvature 
(inflationary) regime the i? 2 -term dominates over the i?-term in the Starobinsky action 
(7), while the coupling constant in front of the i? 2 -action is dimensionless (Sect. 2). 
The Higgs inflation is based on the Lagrangian (160) with the relevant scalar potential 
V A = jA0j (the parameter v is irrelevant for inflation), whose coupling constants £ and 
A are also dimensionless. Therefore, both relevant actions are scale invariant. Inflation 
breaks that symmetry spontaneously. 

The supersymmetric case is similar: the nonminimal action (183) with the X- 
function (189) and the superpotential (188) also have only dimensionless coupling 
constants £ and A, while the same it true for the F (71) -supergravity action with the 
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F-function (190), whose coupling constant / 3 is dimensionless too. Therefore, those 
actions are both scale invariant, while inflation spontaneously breaks that invariance. 

A spontaneous breaking of the scale invariance necessarily leads to a Goldstone 
particle (or dilaton) associated with spontaneously broken dilatations. So, perhaps, 
Starobinsky scalaron (inflaton) may be identified with the Goldstone dilaton! 

The basic field theory model, describing both inflation and the subsequent reheat- 
ing, reads (see eg., eq. (6) in ref. [79]) 

L/V=9 = \d^d^ - V(<f>) + \d^X - \m\x 2 + \lRx 2 + $(ird» - 
- - h(^)<p 

(196) 

with the inflaton scalar field interacting with another scalar field x an d a spinor 
field ip. The nonminimal supergravity theory (175) with the Wess-Zumino superpo- 
tential (187) can be considered as the N — 1 locally supersymetric extension of the 
basic model (196) after rescaling C to (l/\/2)0 c and identifying f = — jj£ because 
of eq. (174). Therefore, pre-heating (ie. the nonperturbative enhancement of parti- 
cle production due to parametric resonance [79]) is a generic feature of supergravity 
models. 

The axion \ an d fermion tp are both requred by supersymmetry, being in the same 
chiral supermultiplet with the inflaton 0. The scalar interactions are 

Mnt(0, X) = ™A0(0 2 + x 2 ) + y(0 2 + X 2 ) 2 (197) 
whereas the Yukawa couplings are given by 

L Y u = \\<KW) + ^AxfeVO (198) 

Supersymmetry implies the unification of couplings since h — — \\ and g 2 = A 2 in 
terms of the single coupling constant A. If supersymmetry is unbroken, the masses of 
0, x an d ^ are a U the same. However, inflation already breaks supersymmetry, so the 
spontaneously broken supersymmetry is appropriate here. 

To conclude, inflationary slow-roll dynamics in Einstein gravity theory with a non- 
minimal scalar-curvature coupling can be equivalent to that in the certain f(R) gravity 
theory. We just extended that correspondence to N = 1 supergravity. The nonminimal 
coupling in supergravity can be rewritten in terms of the standard ('minimal') N — 1 
matter-coupled supergravity, by using their manifestly supersymmetric formulations 
in curved superspace. The equivalence relation between the supergravity theory with 
the nonminimal scalar-curvature coupling and the F(1Z) supergravity during slow-roll 
inflation is, therefore, established. 

The equivalence is expected to hold even after inflation, during initial reheating 
with harmonic oscillations. In the bosonic case the equivalence holds until the inflaton 
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field value is higher than u w M P \/£ hos ~ 10 5 M P \. In the superymmetric case we 
have the same bound u m M P1 /£si4y ~ 10~ 5 M P1 . 

The Higgs inflation and the renormalization group can be used to compute the mass 
of a Higgs particle in the Standard Model by descending from the inflationary scale to 
the electro-weak scale. For example, in the two-loop approximation one finds [80] 

129 GeV <m H < 194 GeV (199) 

with the theoretical uncertainty of about ±2 GeV. It is to be compared to the observed 
bounds known from direct searchers of Higgs particle. For instance, the most recent 
(preliminary) data from the LHC in 201 1 implies 

LHC (ATLAS, CMS) : 124 GeV < m H < 135 GeV (200) 

Therefore, the Higgs inflation is possible in the upper half of the experimentally- 
allowed region (200), though in the lowest portion of the theoretically-allowed region 
(199). It is worth noticing here that in a supersymmetric extension of the SM (like 
the MSSM and NMSSM) there are more particles, when compared to the bosonic SM. 
Hence, the SUSY renormalization group trajectory is going to be steeper, while the 
theoretical SUSY bounds on the Higgs mass at the electro-weak scale are going to be 
lower than those in eq. (199). 



13 Reheating and Quantum Particle Production 

Reheating is a transfer of energy from inflaton to ordinary particles and fields. 

The classical solution (neglecting particle production) near the minimum of the 
inflaton scalar potential reads 

a( * )wfl0 UJ and ^^{sm-J T^To (201) 

A time-dependent classical spacetime background leads to quantum production of par- 
ticles with masses m < oo = M- mi [77]. Actually, the amplitude of ^-oscillations 
decreases much faster [79], namely, as 

exp[~(3H + r)t] (202) 
via inflaton decay and the universe expansion, as the solution to the inflaton equation 

" +3H £ +(m 2 + Tl)ip = (203) 

Here n denotes the polarization operator that effectively describes particle produc- 
tion. Unitarity (optical theorem) requires Im(n) = rriT. The assumption m ^> H has 
also been used here [79]. 
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The Starobinsky model in Jordan frame, 

S = J d^x^fjfsiRj) + S SM (g^ r V) (204) 
after the conformal transformation to Einstein frame reads 

S = ^scalar-tensor gravity^, V 9 ) + S SM (fi^e"^, 1p) (205) 



so that the infiaton ip couples to all non-conformal terms and fields ip, due to the 
universality of gravitational interaction. Therefore, the Starobinsky inflation also has 
the universal mechanism of particle production. 

The perturbative decay rates of the infiaton into a pair of scalars (s) or into a pair 
of spin- 1/2 fermions (/) are given by [21, 22, 84] 



T'l 



M? f M in{ Mj 
r -" = l92% a " d r -"=48^' (206) 

respectively. The perturbative decay rate of the infiaton into a pair of gravitini is [85] 

i ^3/ 2 - 2887T m\ /2 Ml x {ZU,) 

Being proportional to M^ lf , eq. (207) may lead to the cosmologically disastrous grav- 
itino overproduction [86]. However, when the expectation value of the dilaton is of 
the order M P \ (it is the case in the Starobinsky inflation), one can demonstrate that 
eq. (207) reduces to the scalar decay rate (206) proportional to M^ f [86]. 

There is no parametric resonance enhancement here because the produced particles 

rapidly scatter. The energy transfers by the time t rch > (j2 s f T s j^j . One finds the 
reheating temperature [36, 87] 



T - a ]jw^Wr^ m " GeV (208) 

that gives the maximal temperature of the primordial plasma. 

In the context of supergravity coupled to the supersymmetric matter (like MSSM) 
gravitino can be either LSP (= the lightest sparticle) or NLSP (= not LSP). In the 
LSP case (that usually happens with gauge mediation of supersymmetry breaking and 
"^3/2 10 2 GeV) gravitino is stable due to the R-parity conservation. If gravitino 
is NLSP, then it is unstable (it usually happens with gravity- or anomaly- mediation 
of supersymmetry breaking, and m 3 / 2 10 2 GeV). Unstable gravitino can decay 
into LSP. See ref. [88] for a review of mediation of supersymmetry breaking from the 
hidden sector to the visible sector. 
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Stable gravitino may be the dominant part of Cold Dark Matter (CDM) [89]. There 
exist severe Big Bang Nucleosynthesis (BBN) 5 constraints on the overproduction of 
3 He in that case, which give rise to the upper bound on the reheating temperature of 
thermally produced gravitini, T rch < 10 5H ~ 6 GeV [75, 91]. The reheating temperature 
(208) is unrelated to that bound because it corresponds to the much earlier time in the 
history of the Universe. 

When gravitino is NLSP of mass m 3 / 2 ^> 10 2 GeV, the BBN constraints are drasti- 
cally relaxed because the gravitino lifetime becomes much shorter than the BBN time 
(about 1 sec) [75, 91]. In that case the most likely CDM candidate is MSSM neutralino, 
while the reheating temperature may be as high as 10 10 GeV [91]. 

An overproduction of gravitini from inflaton decay and scattering processes should 
be avoided, in order to prevent overclosure of the universe. The cosmological con- 
straints on gravitino abundances were formulated in ref. [86]. Those constraints are 
very model-dependent. 

The rate of decay changes with time, along with the decreasing amplitude of in- 
flaton oscillations. It stops when the decay rate becomes smaller than the production 
rate. Then the particle production accelerates (called pre-heating, or true BB!) due to 
the parametric resonance enhancement [79]. The reheating rapidly transfers most of 
energy to radiation, and leads to a radiation-dominated universe with a oc t 1 ^ 2 . 

In the matter-coupled F(1Z) supergravity with the action 



J d 4 xd 2 9£F(1Z) + R.c. 



after the super- Weyl transformation, £ — > £e , we get 

5* ^scalar— tensor supergravity (-^5 < ^ > ) *^SSM^ (210) 

so that the superscalaron <3> is universally coupled to the SSM matter superfields ^. 



14 Conclusion 

• A manifestly AD, N = 1 supersymmetric extension of f(R) gravity exist, 
it is chiral and is parametrized by a holomorphic function. An F(1Z) supergravity 
is classically equivalent to the scalar-tensor theory of a chiral scalar superfield (with 
certain Kahler potential and superpotential) miminally coupled to the N = 1 Poincare 
supergravity in four spacetime dimensions (with nontrivial G and K), ie. the N = 1 
supersymmetric quintessence. 

The classical equivalence between the F(1Z) supergravity and the quintessence 
N=l supergravity has the same physical contents as the classical equivalence between 
f(R) gravity and scalar-tensor gravity, ie. the same inflaton scalar potential and, there- 
fore, the same inflationary dynamics. However, the physical nature of inflaton in the 

5 See ref. [90] for a review of BBN. 
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f(R) gravity and the scalar-tensor gravity is very different. In the f(R) gravity the 
infiaton field is the spin-0 part of metric, whereas in the scalar-tensor gravity inflaton 
is a matter particle. The inflaton interactions with other matter fields are, therefore, 
different in both theories. It gives rise to the different inflaton decay rates and different 
reheating, ie. implies different physics in the post-inflationary universe. 

Similar remarks apply to the equivalence between Higgs inflation and Starobinsky 
inflation (Sec. 12). The equivalence does not have to be valid after inflation. For 
example, the reheating temperature T reh after the Higgs inflation is about 10 13 GeV 
[73, 74, 76], whereas after the Starobinsky inflation one finds T rch ~ 10 9 GeV [22] — 
see also the next Section. 

By the use of the renormalization group equations in the Standard Model it is 
possible to relate the inflationary scale with the electro-weak scale, and thus 'predict' 
the mass of the Higgs particle. In the two-loop approximation one finds 129 GeV < 
m H < 194 GeV [80]. It is to be compared to the direct searches at LEP: m H > 
114 GeV and LHC: m H < 145 GeV. Therefore, the Higgs inflation is consistent with 
the known experimental bounds the Higs mass in the upper half of the allowed region, 
129 GeV < m H < 145 GeV [80]. In a supersymmetric Higgs inflation, the theoretical 
window for the Higgs mass is lower. 

• It is expected that the classical equivalence is broken in quantum theory 
because the classical equivalence is achieved via a non-trivial field redefinition (Sees. 
3 and 6). When doing that field redefinition in the quantum path integrals defining 
those quantum theories (unter their unitarity bounds), it gives rise to a non-trivial Jaco- 
bian that already implies the quantum inequivalence, even before taking into account 
renormalization. 6 

In the supergravity case, there is one more clear reason for the quantum inequiva- 
lence between the F(1Z) supergravity and the clasically equivalent quintessence super- 
gravity. The Kahler potential of the scalar superfield is described by a. full superspace 
integral and, therefore, it receives quantum corrections that can easily spoil classical 
solutions describing an accelerating universe (those corrections are not under control). 
It was the reason for introduction of flat directions in the Kahler potential and popu- 
lar realizations of inflation in supergravity by the use of a chiral scalar superpotential 
along the flat directions [50, 52, 49]. The F{1Z) supergravity action is truly chiral, so 
that the function F(1Z) is already protected against the quantum perturbative correc- 
tions given by full superspace integrals. It is the important part of physical motivation 
for F(1Z) supergravity. It also explains why we consider F(1Z) supergravity as the vi- 
able and self-consistent alternative to the Kahler flat directions for realizing slow-roll 
inflation in supergravity. 

• The Starobinsky model of chaotic inflation can be embedded into F(1Z) 
supergravity, thus providing the new and truly viable explicit realization of chaotic 
inflation in supergravity, together with the simple solution to the ^-problem, by using 
just supergravity! 

6 See ref. [81] for the first steps of quantization with a higher time derivative. 
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• A simple extension of our inflationary model (Sec. 10) has a positive cos- 
mological constant in the regime of low spacetime curvature (Sec. 11). It is non-trivial 
because the standard supergravity with usual matter can only have a negative or van- 
ishing cosmological constant [82]. It happens because the usual (known) matter does 
not violate the Strong Energy Condition (SEC) [83]. A violation of SEC is required for 
an accelerating universe, and is easily achieved in f(R) gravity due to the fact that the 
quintessence field in f(R) gravity is part of metric (ie. the unusual matter). Similarly, 
the quintessence scalar superfield in F(1Z) supergravity is part of super-vielbein, and 
also gives rise to a violation of SEC. 

In the F(1Z) supergravity model we considered (Sees. 10 and 11), the effective 
f(R) gravity function in the high-curvature regime is essentially given by the Starobin- 

sky function (— ^-R+ ^\ R 2 ). In the low-curvature regime it is essentially given by 

inf 

M 2 

the Einstein-Hilbert function with a cosmological constant, ( ^R — A). Therefore, 

our model has a cosmological solution describing an inflationary universe of the quasi- 
dS type with H(t) = (M 2 n{ /6)(t en d — t) at early times t < t en d, and an accelerating 
universe of the dS-type with H = A at late times. It is the first time when the unifying 
description of inflation and Dark Energy is found in the context of supergravity. 

The dynamical chiral superfield in F(1Z) supergravity may be identified with the 
dilaton-axion chiral superfield in quantum 4D Superstring Theory, when demanding 
the SL(2, Z) symmetry of the effective action. As is well known, String Theory 
supports the higher-derivative gravity. In part, the R 2 A(R) terms may appear in the 
bosonic gravitational effective action after superstring compactification. The apparent 
problem is how to get the anomalously large coefficient in front of the 7£ 3 -term in the 
effective F(1Z) supergravity theory, which would be consistent with string dynamics. 

Supersymmetry in F(R) supergravity is already broken by inflation. The anomaly- 
or gravitationally-mediated supersymmetry breaking may serve as the important ele- 
ment for the new particle phenomenology (beyond the Standard Model) based on the 
matter-coupled F(R) supergravity theory. 

15 Outlook: CP-violation, Baryonic Asymmetry, 
Lepto- and Baryo-genesis, Non-Gaussianity, Tests 

The observed part of our Universe is highly C— and CP— asymmetric (no antimatter). 
Inflation naturally implies a dynamical origin of the baryonic matter predominance due 
to a nonconserved baryon number. The main conditions for the dynamical generation 
of the cosmological baryon asymmetry in early universe were formulated by A.D. 
Sakharovin 1967 [65]: 

1. nonconservation of baryons (cf. SUSY, GUT, EW theory), 

2. C— and CP— symmetry breaking (confirmed experimentally), 



45 



3. deviation from thermal equilibrium in initial hot universe. 

The first condition is clearly necessary. And (in theory) there is no fundamental 
reason for the baryon number conservation. The baryon asymmetry should have orig- 
inated from spontaneous breaking of the CP symmetry that was present at very early 
times, so is the need for the second condition. Then the third condition is required by 
the CPT symmetry, when the CP-violation is compensated by the T-violation, so it 
has to be no thermal equilibrium. 

There exist many scenarios of baryogenesis (see ref. [66] for a review), all designed 
to explain the observed asymmetry (BBN,CMB): 

(3 = nB ~ n B = ( 6 .o ± 0.5) • 1(T 10 (211) 

Here ub stands for the concentration of baryons, n-g for the concentration of anti- 
baryons, and n 7 for the concentration of photons. 

Perhaps, the most popular scenario is the nonthermal baryo-through-lepto-genesis 
[67, 68], ie. a creation of lepton asymmetry by L-nonconserving decays of a heavy 
(m w 10 10 GeV) Majorana neutrino, and a subsequent transformation of the lepton 
asymmetry into baryonic asymmetry by CP-symmetric, B-nonconserving and (B-L)- 
conserving electro-weak processes. 

The thermal leptogenesis requires high reheating temperature, T rch > 10 9 GeV 
[92], which is consistent with eq. (208). 

The matter-coupled F(1Z) supergravity theory may also contribute towards the ori- 
gin and the mechanism of CP- violation and baryon asymmetry, because 

• complex coefficients of F(JZ) -function and the complex nature of the F(1Z) 
supergravity are the simple source of explicit CP- violation and complex Yukawa cou- 
plings; 

• the nonthermal leptogenesis is possible via decay of heavy sterile neutri- 
nos (FY-mechanism) universally produced by (super)scalaron decays, or via neutrino 
oscillations in early universe [93]; 

• the existence of the natural Cold Dark Matter candidates (gravitino, axion, 
inflatino or, maybe, inflaton itself!) in F(1Z) supergravity; 

• F (1Z) -supergravity can naturally support hybrid (or two-field) inflationary 
models because it already has a pseudo-scalar superpartner (axion) of inflaton. As 
is well known, non-Gaussianity is a measure of inflaton interactions described by its 
3-point functions and higher - cf. eq. (24). The non-Gaussianity parameter / NL is 
defined in terms of the (gauge-invariant) comoving curvature perturbations as 

n = n gI + ^/nl^p (212) 

The non-Gaussianity was not observed yet, though it is expected. As regards the 
single-field inflationary models, they predict [94] 

/NL = 4(l-n a )«0.02 (213) 
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Finally, we would like to comment on possible testing of f(R) gravity and F(JZ) 
supergravity in Solar system and ground-based experiments. 

As regards the large-scale structure of the present universe, the scalaron (ie. the dy- 
namical spin-0 part of metric) may be responsible for its acceleration or Dark Energy 
(Sec. 11). However, since scalaron is universally coupled to all matter with gravi- 
tational strength, it may easily lead to a large violation of the equivalence principle. 
To avoid that, the scalaron should be essentially massless on the Solar system scales, 
because of the strong observational constraints from experimental tests of the equiva- 
lence principle [95, 96]. Moreover, it also should not give rise to a large violation of 
the equivalence principle in ground-based (on Earth) laboratories, because of the tight 
constraints on the fifth fundamental force in Nature [97]. 

A natural solution to both problems is provided by Chameleon Cosmology [98, 
99], where the effective scalaron mass is dependent upon a local density p (see also 
refs. [100, 101]). The scalar potential of the scalaron (Chameleon) field gets an extra 
term, which gives rise to an effective scalar potential of the form 

V eS (<p) = V(ip) + pexp (P<p/M Pl ) (214) 

where the parameter (5 is of the order 1. The exponential factor here arises due to 
the universal coupling of the scalaron to other matter of density p — see eq. (205). 
Therefore, when the effective mass is suffuciently large in a dense environment, one 
can evade the known constraints on the equivalnce principle and the fifth force, since 
the resulting violations are exponentially suppressed. 
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